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We describe packages for the calculation of radiative corrections to two-fermion production at colliders. The 
packages use DIANA, and also QGRAF, FORM, Fortran, and further unix/linux tools. The one-loop calculations in 
the Standard Model are highly automatized with the package aiTALC. Further, the automatic determination of 
all the matrix elements for two-loop corrections to massive Bhabha scattering in QED and the classification of 
their topologies and prototypes is done with DIANA. A generalization to the Standard Model is straightforward. 



1. INTRODUCTION 

For the future e+e~ Linear Collider (LC) we 
need quite precise predictions of cross-sections for 
a variety of reactions, among them 
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The predictions have to be calculated with ac- 
count of quantum corrections, typically with one- 
or two-loop accuracy. This has to be done in some 
model of field theory, notably the Standard Model 
or the Minimal Supersymmetric Standard Model. 
First complete one-loop calculations in the Stan- 
dard Model for two-fermion production pQ and for 
Bhabha scattering [2] date back to 1979. After 25 
years of perturbative calculations, the quest for a 
more or less complete automatization of this kind 
of calculation is natural. Nevertheless, only few 
packages for such a task are publicly available. 
We would like to mention CompHEP [3I4I5I6I7] . 
FeynCalc/FeynArts/FormCalc/LoopTools/FF 
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|8I9I10I11| . the grace project |12I13| and SANC 
|14| . These packages have a variety of options to 
be used, but they do not go substantially beyond 
the one-loop level. 

We are interested in a package for one- and two- 
loop calculations and we think that a prospective 
approach might be based on DIANA )15I16I17) . an 
interface to QGRAF 18;,, to be combined also with 
FORM |19I20| and Fortran, c++ and additional 
unix/linux tools. 

Recently, we performed complete, high- 
precision one-loop calculations for (/ ^ e) 
in the Standard Model |21I22| and demonstrated 
numerical agreements with the results of other 
groups, finally with up to eleven digits |23| . In 
a next step, the program package was completely 
rewritten in order to allow also the unexperienced 
user to create their own numerical code. The re- 
sult is the package aiTALC. It covers also Bhabha 
scattering and will be described in more detail in 
the next section. In an earlier comparison for 
Bhabha scattering |24j a per mill agreement was 
established (see also 25 2TJ). For realistic ap- 
plications, one has to include also higher order 
corrections and to combine the so-called 'weak li- 
brary' with a Monte Carlo code for the treatment 
of real bremsstrahlung; this is not discussed here. 
A dedicated comparison of this kind proved an 
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accuracy of the order of 10 4 |27j . 

2. DIANA and two-fermion production 

It is not so long ago that DIANA [2Hj was in 
a kind of experimental state. There exist sev- 
eral variants of the package with identical version 
numbers but with quite some different properties. 
For this reason, we decided to install a chain of 
versions with well-defined version properties at 
our location, presently the last one being v. 2. 35 
|29|. From the user's point of view, the package 
is characterized by a small number of input files: 

• process. cnf - the file defines how to run 
DIANA; for examples see [22]; the incom- 
ing and outgoing particles plus loop order 
should be edited; 

• Model . model - we have prepared four model 
files: QED. model, StandardModelO .model 

(a basic file with leptons 

only), StandardModel .model (com- 
plete with all masses and mixings), 
StandardModelC. model (with counter 
terms for neutral current 2/ — > 2/, fol- 
lowing PHI, a typo corrected with [3"T]1: 

• FeynmanRules . f rm we have pre- 
pared three sets of Feynman rules: 
FeynmanRulesO . f rm (a simple version), 
FeynmanRulesl . f rm (with account of 
actual external momenta distribution), 
FeynmanRules2 . f rm (with account of in- 
ternal momenta). 

These files, together with a proper application 
of QGRAF and of additional DIANA options, give 
a huge flexibility to prepare a sample of prepared 
matrix elements of a given loop order in a given 
model for a subsequent FORM or MAPLE calcula- 
tion. As an application in the Standard Model, 
the flavor non-diagonal reactions 
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may also be treated. In this respect it would be 
extremely useful to have additional model files, 
e.g. for the minimal supersymmetric Standard 
Model. Some time ago, a project in XML was 



proposed, being intended on creating a standard 
for model files to be used in different CAS |32|. 

2.1. Full automatization of one-loop 
EWRC: alTALC 

A completely automated example, based on 
DIANA and other packages (see Fig. QJ, is the 
project alTALC [HSj) that provides also numeri- 
cal output. 



alTALC 




Figure 1. Flow chart of the alTALC package. 



It is designed to organize the calculation of 
tree-level and one-loop corrections to 2 — > 2 
fermion processes, including Bhabha scattering 
(and in the near future also Moller scattering). 
This tool allows presently the calculation of unpo- 
larized differential and integrated cross-sections. 
Differing from topf it |34|. alTALC treats all par- 
ticle masses (including the electron mass) and 
mixings exactly, but does not cover hard photon 
emission. 

In a simple driver file process . ini (simple 
form of process . cnf), the user specifies the ingo- 
ing and outgoing fermions. The determination of 
matrix elements, form factors and other kinemat- 
ical routines proceeds then automatically. Later 
on, two Fortran files give access to the specific 
values of the parameters in our model and the 
numerical output of the code (i.e. a number of 
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Table 1 

Cross-sections for e~e + — > e~e+ (7) at y/s 

500 GeV with a photon energy cut i?" la f * 
V^/H). 



4. MASSIVE TWO-LOOP 
SCATTERING 



BHABHA 



COs6> 



d cos # 
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-0.9 0.19344 50785 26863 6 aiTALC 

-0.9 0.19344 50785 26862 2 FA/FC/LT 

-0.9 0.19344 50785 62637 9 m e = 

0.0 0.54667 71794 69423 1 aiTALC 

0.0 0.54667 71794 69421 8 FA/FC/LT 

0.0 0.54667 71794 99961 4 m e = 

0.9 0.17292 83490 66507 2 • 10 3 aiTALC 

0.9 0.17292 83490 66508 • 10 3 FA/FC/LT 

0.9 0.17292 83490 61347 4 • 10 3 m P = 



data points in the angular distribution, integrated 
cross section, running flags, etc). 
Advanced features are not extremely user 
friendly, but still under development. Fermion 
masses are included by default, but may be ne- 
glected, and soft photonic corrections may be in- 
cluded (or not), as well as one may perform a full 
one-loop tensor integral reduction to the master 
integrals Aq, Bq, Co and Dq in the Passarino- 
Vcltman scheme £Q. All these possibilities will 
be described in more detail in a tutorial to be 
published soon. 

Some more details may be found also in |23| . 

3. NUMERICAL ONE-LOOP RESULTS 
FOR BHABHA SCATTERING 

The package aiTALC was used for a pre- 
cise calculation of the one-loop electroweak 
corrections to fermion pair production and 
Bhabha scattering at LC energies. We com- 
pared the result of this with numerics from 
FeynArts/FormCalc/LoopTools 

The Table Q] shows an agreement of the two 
calculations of 14 digits for the 0(a) corrections 
(i.e. Born+QED+weak+soft) for the case with 
an exact treatment of m e ^ 0, while the simplified 
comparison with m e — * does not provide more 
than 10 digits agreement. For practical purposes, 
this make no difference, of course. 



In a related activity, we use DIANA for a prepa- 
ration of two-loop matrix elements in massive 
QED. We apply it for the solution of two different 
problems: 

• The calculation of interferences of the two- 
loop amplitudes with the Born amplitudes 
(i.e. of the integrands for the determination 
of the scalar integrals); 

• The determination of all prototypes. 

Prototypes are topologies of diagrams, where ad- 
ditionally the different masses on internal lines 
are taken into account. An example is shown in 
Figure[3 where all prototypes with three external 
lines and six internal lines for the Bhabha pro- 
cess are given. Diagrams V1-V5 are genuine two 
loop QED vertices for the Bhabha process. Ver- 
tices V6-V10 come from extractions of one sin- 
gle line from two-loop QED box prototypes (al- 
together there are six genuine box prototypes). 
Strictly speaking, the prototypes V4-V6 and V8 
should already be shrinked to five internal lines 
and classified into this class of prototypes because 
they have two internal lines with identical mo- 
mentum. There are many more prototypes with 
five internal lines which come from the extraction 
of one (or two) lines of original vertex (or box) 
prototypes; we do not present them here. For 
these cases and for cases where even more lines 
are extracted, an automatization of the procedure 
for the identification of prototypes is very useful. 
Technically it has been made by merging informa- 
tion from DIANA and QGRAF in a subsequent FORM 
program. 

Let us mention finally that out of all the proto- 
types shown in Figure |3 only V2 and V10 corre- 
spond to master integrals; all the others may be 
reduced to other master integrals. For a complete 
set of master integrals see ■ 
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Figure 2. Vertices for the two- loop Bhabha pro- 
cess with six internal lines. 
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